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Abstract 

We apply the relativistic Bethe-Salpeter formalism to the calculation of final- 
state effects in the production of a tt pair at threshold. We find that final-state 
rescattering does not affect the momentum distribution of the tt pair to lowest 
order in the strong coupling constant. This result correctly extends earlier results 
based on the non-relativistic Lippmann-Schwinger equation. 
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As was first pointed out by Fadin and Khoze |]^, the structure of the tt threshold region is 
very interesting. The cross-section is enhanced by tt resonances and depends strongly on the 
mass of the top quark m*, its decay width F^, and the strong coupling constant as{mt). Unlike 
the other heavy quarks (b and c), tt pairs cannot form narrow resonances. This is due to the 
large value of rrit, which implies that the dominant decay mode is the weak decay t W^b 
with a width F^ ~ 1 GeV. Thus F^ exceeds Aqcd and as a result the tt pair decays before it has 
time to hadronize, which makes perturbation theory applicable. Since theoretical predictions 
are possible, it would be extremely interesting to measure the cross-section for threshold tt 
production at the next linear e’''e“ collider (NLC). 

In the threshold region, the tt pair is produced with little kinetic energy. Therefore, one can 
employ non-relativistic techniques to calculate the cross-section. In analogy with the Hydrogen 
atom, the tt pair is (weakly) bound together by a QCD Coulomb-like effective potential. The 
cross-section for tt production is related to the Green function, which obeys an appropriate 
Lippmann-Schwinger equation, by the optical theorem |^. Although it is possible to use this 
method to calculate higher-order corrections to the cross-section [^, there are subtleties which 
ought to be taken into account. As was shown by Kummer and Modritsch [^, these problems 
can be avoided if one employs the relativistic Bethe-Salpeter formalism, in analogy with the 
abelian positronium case 0]. 

In ref. §, we extended the results of ref. 0 by showing that electroweak corrections to 
the decay width of toponium are suppressed by at least four powers of the strong coupling 
constant, To this end, we employed the covariant Lorentz gauge and perturbed around the 
solution to the Bethe-Salpeter equation in the instantaneous approximation. The calculation 
was manifestly gauge-invariant and correctly took into account all contributions to the decay 
rate (Goulomb enhancement and phase-space reduction effects). 

Here, we wish to apply the argument of ref. to the calculation of radiative corrections 
to the momentum distribution of the tt pair produced at threshold at a future e’''e“ collider 
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(e.g., the NLC). One-loop corrections to the total cross-section of tt threshold prodnction have 
already been shown to vanish to hrst-order in by Melnikov and Yakovlev 0 , who employed 
the non-relativistic potential formalism. By using the relativistic Bethe-Salpeter formalism, we 
shall show that there are no o{as) corrections to the differential cross-section either. 

We are interested in one-loop corrections to the process (Fig. P 

e+e"—^7,^°— >ti — >W+bW-b, (1) 


where at threshold the dominant contribution is due to the tt bound states (toponium). To 
simplify the discussion, we shall ignore tt production via the boson. We shall also neglect 
the Higgs boson exchange between t and t, which contributes to the formation of bound states. 
Our calculation can easily be extended to accommodate these additional effects. 

Let us first review the non-relativistic potential formalism in order to demonstrate its power 
as well as its limitations. The tt system at threshold forms short-lived bound states which are 
eigenstates of the Hamiltonian 


M, 


Cptts 


+ 2Mt , 


where we have included the effects of the decay in the mass parameter Mf 
Green function for this Hamiltonian, 


( 2 ) 

rut + iVt- The 




fe(x)Vi;(x') 
E - En- iT„ 


is related to the cross-section of tt production by the optical theorem. 


a(7 ^ tt) = 


mf 


Im G(0, 0; E) , 


( 3 ) 

( 4 ) 


where E = s/s — 2mt. This expression for the cross-section includes the contributions of all 
ladder diagrams containing an arbitrary number of Coulomb-like gluon exchanges between the 
two quarks (t and t). This class of diagrams dominates in the threshold region. By splitting 
the Hamiltonian (H) in its real and imaginary parts, El = Hq + iT, and expanding. 
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(5) 


:iT- 


H - E Ho-E Ho-E Ho - E 
we may write to first-order in the Fermi constant [| 

d?p 


-I- . . . , 


a(7 -^tt) = 




ml 


(2vr) 


\G{p-,E)\^T{p-,E) . 


( 6 ) 


The fnnction r(p; E) represents the decay rate of toponinm. To lowest order, T = 2Tt. Higher- 
order corrections are dne to time dilatation and are negligible [|,i- Thus, the momentum 
distribution of the top quark to lowest order is given by 

da “iQlo? 


d\p\ 


2mf 


V |G(p;i?)h 


(7) 


Higher-order corrections to the cross-section (Q) can be calculated in this formalism using 
perturbation theory. However, one needs to be careful in applying ordinary perturbation theory, 
because the decay products (e.g., b and b) feel the binding QCD potential and therefore, they 
cannot be described by plane waves. To avoid such subtleties, one may use the relativistic 
Bethe-Salpeter formalism which is manifestly gauge-invariant. 

To calculate the differential cross-section, we note that it factorizes into two parts, describing 
tt production ^ 7 , ^ ti), and scattering and decay of the tt pair (ti bW~^ bW~), 

respectively. We shall concentrate on the latter part. Let be the total incoming momentum, 
and be the relative momentum of the tt pair. After we integrate over the hnal-state phase 
space, we can use the optical theorem to express the cross-section as 

da 


d\p\ 


ImAti{k,p,p') , 


( 8 ) 


where Au is the amplitude for tt scattering, and p^ is the relative momentum of the hnal tt 
pair. This four-point amplitude satishes the inhomogeneous Bethe-Salpeter equation. 


iThe absorptive part of the Hamiltonian, F, is proportional to Gp. 
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( 9 ) 


Y[^^\p+)Yi^‘^\p_)Ati{k,p,p) = 1 + / ^^V{p,p'';k)Ati{k,p',p') , 
where 11 (p) is the complete inverse fermion propagator, 

n(p) =^-- E(p) , (10) 

and we have defined momenta 

P± = ^±P- (11) 

V{p, p'] k) is a potential function which consists of the two-fermion irreducible graphs. For our 
purposes, the mass is complex. 


Mt = mt + iVt . 


( 12 ) 


The solution to Eq. (||) has poles which are due to bound states. Near a pole, 

kXk{p)Xk{p') 


Att{k,p,p') 


k‘^-M^ 


(13) 


The bound-state wavefunctions Xfc(p) satisfy the homogeneous Bethe-Salpeter equation, 

U^^\p+)U^^\p_)Xkip) + j -^^V{p,p';k)xk{p') = 0 . (14) 

which is represented graphically in Fig. p. 

To lowest order in ag and neglecting electroweak interactions, the potential is 


Va(p, p']k)^CF 410.7^(3'“'(p - p')7<^> , 


(15) 


where Cp = 4/3 is the Casimir operator, and G^'^{q) is the lowest-order gluon propagator. In 
the Feynman gauge (omitting group theory factors). 




9 I ' ’ 

+ ^e 


( 16 ) 
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and the potential Vo{p,p'-,k) is independent of k^. At threshold, the quarks move with 
non-relativistic velocities and the Bethe-Salpeter equation can be approximated by the non- 
relativistic Schrodinger equation in momentum space, and then solved. To this end, we shall 
work in the total rest frame in which the overall momentum is k^ = [E, 0). In the instantaneous 
approximation, the potential becomes 


Vr\p,p';k)=CFA7rasil, 


( 1 ) 


r76 


( 2 ) 


(17) 


{p-p'y 

If we integrate over we can write the Bethe-Salpeter equation (|^ in terms of the wave- 
function <h(p) = / ^xip) 

d?p 1 


(//•'> + - At’A?*)CF4ira, y 


. ( 18 ) 


(27r)3 (p — p Y 

where H is the Dirac Hamiltonian and A_|_ (A_) is the projection operator onto positive (neg¬ 
ative) energy states. In the non-relativistic limit, this reduces to the Schrodinger equation in 
momentum space 


Mt 


+ 2Mt — Ej <h(p) = C® dvrcTs J 

Thus, we obtain the energy levels 

MtCla^ 


(Ep 


(2vr)3 {p-p ) 


t^Mp) 


(19) 


En, — 2 Mi 




+ o(aJ , 


( 20 ) 


which are the Bohr levels of the Coulomb-like QCD potential (ITtI). Therefore, the hrst-order 
QCD correction to the decay rate of toponium is 

^ OWE 


r*+ — 2r+ 1 — 


8n? 


( 21 ) 


which may be attributed to time dilatation The spherically symmetric S' = 0 states are 
given by 




,(p) = (M.CFa.)-^/ \YEY 


(1 -I- n2|/)"-+i 


y = 


4:P‘ 


MfClal ' 


( 22 ) 
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where is a polynomial of order n — 1 related to the Laguerre polynomials. For n = 1, we 
have Cl = 16a/^. 

To discnss hnal-state corrections, we need to consider the diagrams in Fig. |^. They give 
hrst-order corrections to the differential cross-section for tt decay, according to the nnitarity 
theorem. This may be seen by cntting these diagrams to prodnce the graphs of Fig. ^ In 
the positroninm case, these loop graphs represent magnetic moment effects and contribnte to a 
o{a^) shift in the energy levels (poles). We have shown that there is no W boson contribntion 
to the color magnetic moment of the top qnark p . This is due to the fact that W only couples 
to a left-handed current. We shall sketch the proof of this statement for completeness. 

The energy level shifts can be calculated by perturbing around the solution to the 


Schrodinger equation (|T^. The potential to be treated perturbatively is V — Vf- 


inst 


There 


is also a contribution from the disconnected diagrams which are due to the self-energy terms in 
the fermion propagators (Eq. ([T0|)), but they can be absorbed in the potential if we make use of 
the Schrodinger equation. Thus, according to the Bethe-Salpeter formalism [Q, the hrst-order 
energy level shift is 

AEn = ($n| D,ip) _ E^) {V - vr') - E^) D,{p) |<F„) , (23) 


where the inner product involves an integral over the four-momentum. H is the Dirac Hamil¬ 
tonian, and Dk is the product of two free fermion propagators (c/. Eq. (pH))), which can be 
expressed in terms of the projection operators A± as 

aL"aL"> 


De(p) £ +p“ ± (Ep - m)J|£;/ 2-pO ± (Ep - ie)] ’ 


(24) 


where Ep = \Jp‘^ + rril is the energy of the quark on the mass shell. 

To lowest order, the potential is Vq — (Eqs. ( p!5|) and ([T7|)). This is analogous to the 
positroninm case, and produces an o(a^) shift in the energy levels. The hrst-order electroweak 
correction is 
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Vi{p,p'; k) = 47rCi.a,aty “ p')^y\p-^p'-)^ , (25) 


where p± = k /2 ± p, p'^ = k /2 ± p' , and we have made explicit the electroweak coupling 
constant aw ~ GpM^, where Gp is the Fermi constant and Mw is the mass of the W boson. 
The vertex function A^{p,p') consists of the diagrams shown in £g. |^. It is guaranteed to give 
a gauge invariant contribution by the Ward identity satished by the one-particle irreducible 
function, 

{p-p'n,{p,p') = n{p)-U{p’) . (26) 

Since we are only interested in hrst-order corrections, we may replace by its real part 
The contribution of Vi to the energy level shift (Eq. (^3])) can then be written as 

^ C\alaw r d^p d^p' _ Gn{n^y) Gn{n^y') 

”■ IQrrit J {2n)^ (271)^ {p — p'Y + ie {1 + n'^y)^ (1 + n'^y')^ 

X (^Dkip')(^A\l\p+,p'^hl^'^+ -f'd^Al^\p-,p'_)jDk{p)'^ , ( 27 ) 

where y = Ap'^/mfCpal and y' = Ap'"^/mfCpal- A simple scaling argument shows that the 
lowest-order contribution to the integral comes from the small three-momentum region. Momen¬ 
tum insertions contribute additional powers of as- At low momentum transfer, the three-point 
vertex A^ may be written in general as 

Mp,p') = + (T^.uk'^d^2{k^) , ( 28 ) 

where k = p — p', and = |[ 7 ^, 7 ^]- The form factors and J^2 are regular as k"^ —>• 0. 
In the positronium case, T2 gives an o(a^) contribution to the energy level shift, and is due to 
the magnetic moment interaction. In our case, we need to multiply the gamma matrices by the 
projection operator |(1 — 75 ), due to parity violation of weak interactions. A straightforward 
explicit calculation shows that the form factor J^2{k‘^) vanishes to lowest order in k‘^. It follows 
that the three-point vertex is proportional to {p — p'Y (recall that = P- — = p — p')- 
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Having established the leading-order behavior of A^, we can now estimate the integral in 
Eq. (1^). As we jnst showed, A^ contributes a factor {p — p'Y- This factor cancels the gluon 
propagator. Then the integral over pq and pg can be easily done, because of the respective poles 
in the operators D{p) and D{p'). The resulting expression contains six three-momentum factors 
implying that the integral is o(a®). Therefore, the electroweak correction to the decay width is 
negligible. Of course, no conclusion can be drawn regarding the exact value of the electroweak 
correction, because such a high order is beyond the scope of hrst-order perturbation theory. 

Having shown that the poles do not get shifted due to hrst-order hnal-state corrections, 
we deduce that there are no hrst-order corrections to the amplitude either. Therefore, the 
total cross-section as well as the diherential cross-section of tt production and decay do not get 
corrected to hrst-order by hnal-state rescattering. 

In conclusion, we have presented a simple physical argument showing that hnal-state thresh¬ 
old ehects vanish to hrst-order in perturbation theory. This is true for the diherential cross- 
section of tf production, which strengthens a previous result on the total cross-section by 
Melnikov and Yakovlev |^. Our argument was based on the relativistic Bethe-Salpeter for¬ 
malism. No special gauge-hxing procedure was required and the calculation was manifestly 
gauge-invariant. 
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FIGURES 




FIG. 2. The homogeneous Bethe-Salpeter equation. 
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FIG. 3. First-order electroweak corrections to tt production. 




11 














